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Abstract 

We provide a correction to the Stefan-Boltzmann law and discuss the problem 
of a phase transition from the superfluid state into the normal state. 
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In his book [I], Schrodinger notes that the asymptotic form at a large number of 
particles N can be obtained with a precision not better than \/N . The leading term of 
the asymptotic form as N — > oo with the volume V — > oo but with N/V — > const is 
called the thermodynamic limit. In the example of the Stefan-Boltzmann radiation law, 
we show that the next asymptotic term has the order N 2 ^ 3 and that it is impossible to 
achieve a better precision. We also give this term explicitly. We say that this term, i.e., 
the limit of the difference between the exact answer and the thermodynamic limit, divided 
by iV 2//3 , is the quasithermodynamic limit. 

In this note, we present some general considerations about the thermodynamical limit 
in statistical physics. First we dwell on the study of the notion of black body. 

The Rayleigh- Jeans formula describing black bodies in classical physics, valid for low 
frequencies, was extended in 1900 to high frequencies in the form of the famous formula 
due to Planck, who proposed to consider discrete energies and introduced the constant 
that now bears his name. That formula implies, in particular, the Stefan-Boltzmann 
formula, which had been discovered earlier. 

The derivation of the Rayleigh- Jeans formula is based on the Maxwell equation and 
the Gibbs distribution in classical mechanics. It also makes use of the complete isotropy 
of black emission. 

On the basis of the assumption that black body emission is completely isotropic (see 
[2], p.205]), we rigorously obtain a correction to the Planck formula and to the Stefan- 
Boltzmann law. In many precise experiments about black radiation, the necessity of 
this correction was not noticed. On the other hand, the exact Maxwell equation for free 
photons, which is used to derive the Rayleigh- Jeans law, does not contain the parameter 
h, and the appearance of a correction to a law containing this parameter, leads to the 
problem of carrying over the correction containing the parameter to the Maxwell equation. 

As to the Stefan-Boltzmann law, which reads 

F = -^VT\ (1) 

3c 

where F is the free energy, V is the volume, c the speed of light, and 5 the Stefan- 
Boltzmann constant, 5 = 6 q Jl c s , it is a fact that here small discrepancies between theory 
and experiment were observed. Thus, the Physical Encyclopedic Dictionary of 1966 [3], 



1 



p. 82 says: "The experimental value of 5 exp is somewhat larger than its theoretical value 
5theor obtained by integrating the Planck formula over the wavelength A (or the frequency 
v). The reasons for this discrepancy are not quite clear." 

The quasithermodynamic correction to this law is given by 

F = -|W-2C(3)^V 2 / 3 , (2) 

where ( is the Riemann zeta-function. This result follows from Theorem 3 in pertaining 
to number theory. This general theorem is applicable to objects in economics, linguistics, 
and semiotics. In our case, it determines the domain where most of the choices are 
positioned that satisfy the inequality 

^2(i + j + k)N iJtk <E, (3) 

where E ^> 1 is fixed, and the Nij^ are any natural numbers (the number of choices 
corresponding to ([3]) is finite). The theorem is roughly formulated as follows (see [5] for 
the exact formulation): the probability that the number of choices - 

i yj {i+j+mj:t - £ j+i+^ >ENi», ( 4) 

i+j+k<s i+j+k<s 

where s is fixed, is exponentially small as E — > oo; in other words, there are not too 
many choices outside the specified interval. The parameter b is here determined from the 
equation 

2.^1 e (i+j+k)b _ l - C/ ' W 

and N is the mean value of the sum £\ . k N iy j t k under condition ($\), E = E/N. 

This theorem can be applied to quantum statistics if the following assumptions are 
made: 

1 . All choices of particle distribution over energy levels are equiprobable (something 
of the sort of the equidistribution law, which we call the "no-preference law"). 

2. Blackbody radiation is totally isotropic. 
Suppose we have the condition 

J2(i+j + k)N ijj>k <E, (6) 

where E is a constant, Nij^ are "equally chosen" integers, i.e., are equiprobable or equally 
distributed. Then it follows from the analog of the theorem in [1] that most of the variants 
will cluster near the following dependance of the "cumulative probability" : 



l 

Bi= Y N hi,k, 

i+j+k<l 



with I < s, s = max(l + j + k) in inequality 



Bl ~2^Y e «-l ' (7) 

i=i 



where b is determined from the conditions 

^ 2 e u - 1 v ; 

if — >■ oo. Denote by N the mean value of the sum £V . k Nij jk under condition (jSJ). 

It turns out that, in this case, most of the variants will also cluster around the following 
dependance of "local energy": E\ = ^2 i+ j +k< i(i + j + k)N it j ik , E[ < E, near 



1 z(i + l)(i + 2) 



,6/ 



(9) 



8=1 



where b is determined from conditions ([5]) if £7 — > oo. 

The theorem from [5] can then be stated similarly to the theorem from [I]. 

Theorem 1 Suppose that all choices of the families {iVj} such that 

^2(i + j + k)Nij tk <E, (10) 

are equiprobable. Then the number of variants Af of families {Nij jk } satisfying (0|) and 
(Td\) as well as the following additional condition: 



(i+j + k)N i>jjk -Y,? 



1 £(i + l)(i + 2) 



obi 



> EN 1 ' 2 (In iV) 1/2H 



(11) 



i+j+k<l 8=1 

is less than {c\M)/M^, where c\ and m are arbitrary numbers); here E — E/N ~ 1/b. 

Notation: M. is the set of all families {Nij tk } satisfying condition (jH]), Af{A4} is the 
number of elements of Ai, A is a subset of A4 satisfying the condition 



^1 i(i + l)(i + 2) 



i=0 



1 



< A, 



I = 0,1,..., s, 



where A and b are some real numbers not depending on /. 
Denote 



1 i(i + l)(i + 2) 



i=i 



Let us recall the scheme of the proof, similar to that of Theorem 3 from in our case. 
It is obvious that if Ai is the number of families {N{j jk }, then 

M{M \ A} = 

s 

E { @ i E " E( i+ ^' + ^ v ^}^>: ! ^-,^'X i , .vn B {^ - A l})- ( 12 ) 



{ N i,j,k} 



1=0 



Here the sum is taken over all integers Nij jk , 0(A) is the Heaviside function, and 5 kuk2 is 
the Kronecker delta. 

Using the integral representations 
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8nn> = — / dipe 

Z7T 



-iN<p iN'tp 



by(l+i\) 



i r°° i 
e(y) = — d\ — 

2m J.oo A - i 
and performing the standard regularization, we obtain 

dEQ^E- £(i + J + k)N hhk ) e 



i,j,k 



e -6Ei,i,fc(*+i+fe)-^i,j,* 



Denote 



Z(6,JV) = e_&El 



j,fc( i +J+ fe ) JV i,J,fe 



where the sum is taken over all Nij^. Further, introduce the notation 



' 1 
Ci(ia,b) = JJ&(ia,&), &(iot, 6) = ( _ ,_ b 

8=1 ^ 



bjb'0'+l)/2 ' 



i = l,...,Z. 



It follows from (|T3|) that 



hence 



Z(b,N) = ±- £ dae- lNa ( s (ia,b); 



M{M \A}< 



z(2?r) 



X 



(13) 
(14) 



(15) 



(16) 



exp(— iA^y?) 

xIJed^-Ai) 



(17) 



where b is a real parameter for which the series converges. 

Estimating the right-hand side, carrying the absolute value sign under the integral 
and then further under the sum, we obtain after integration over ip 



bE 



M{M\A}<— exp{-by2(i + j + k)N i>jtk } 

{Ni,j,k} ij,k 



xn©d^-Ai). 

1=0 

From the inequality for the hyperbolic cosine cosh(x) = (e x + e~ x )/2 

s s 

2 s JJcosh(xj) > e s Vx/: ^ \x t \ > 5 > 0, 



(19) 



it follows that, for all positive c and A, we can write (compare [6], [7]) 



IJeflSi- A|) <Te- cA f[cosh(c ^ (i + j + k)N i>j;k - c0 6 Y 



(20) 



1=0 1=0 v i+j+k<l 

where 

^1 i(i + l)(i + 2) 



2 e bi - 1 

i=l 



Thus, we obtain 



(21) 



M{M \A} < e- cA exp{bE) ^ exp j-6^(i + j + k)N ijj)k \ 

x Y\ cosh ( ^ c(i + j + k)N id)k - ap b J 

Z=0 ^i+j+k<l ' 

s 

= e bE e~ cA n(0(0, & - c) exp(-c^) + 0(0, b + c) exp(c^)}. 

1=0 

Let us apply Taylor's formula to 0(0, b ± c). There exists a 7 < 1 such that 
ln(0(0, 6 ± c)) = In CM 6) ± c(ln0)U0, 6) + y (hi0)l'(0, 6 ± 7 c). 

Obviously, 

|ln0--^. 

Let us put c = A/D(0, b), where -D(0, 6) = (In 0)&(0, b) is positive for all b and monoton- 
ically decreases as b increases. The right-hand side of (12"T1) does not exceed 

As in [4], we obtain 

AT(M \ A) < e bE d(0, b)e~ £A2/D(0 ' b) . (22) 

Therefore, in the interval from — 7r to 71 over which the integral (|T6|) is taken, the only 
contribution comes from a neighborhood of the point a = 0. Let us compute the integral 
( |T6l) by the Laplace method with precision up to N~ m and then apply all the subsequent 
arguments from the proof of Theorem 3 in [1] . 

In order to estimate 0(0, (3) from below, we can use the exact asymptotics obtained by 
Krutkov [8] for the three-dimensional oscillator. In the case of a general spectrum A n , the 
formula obtained by the saddle-point method meets with considerable difficulties due to 
the fact that an infinite number of saddle points can appear in certain concrete examples 
(Koval', private communication). 

Note that the sum over all variants satisfying inequality ([6]) may be interpreted as a 
discrete continual "path integral", the paths being the variants. The asymptotic leading 
term of the continual integral is concentrated, as a rule, near one principal "trajectory" 



(this is the Laplace method for the continual integral), and rapidly tends to zero outside 
a neighborhood of this trajectory. This is expressed by the statement of the theorem. 
The size of this neighborhood (in our case EO(N 1 ^ 2 (lnN) 1 ^ 2+e )) determines the limiting 
precision with which it makes sense to compute this "principal trajectory" . 

In this connection, we propose the following somewhat modified conjecture about the 
Schrodinger rule that he qualified as a "law of nature" . If a certain number of particles, 
molecules, genes in a chromosome is equal to N, then one can obtain a statistical law with 
precision of no more than 0(y/N In N). The mathematical meaning of this conjecture is 
that estimate (Hip cannot be improved by more than e. 

Since E — > oo in ([S]), it follows that we can put s = oo in the sum 

^ 2(e bi - 1) ' 
i=i v ' 

we can apply the Euler formula of the form 



E /(")= / f(x)dx + p(b)f(b)-p(a)f(a) 

b>n>a ^ a 



+ cr(a)f'(a) — a (&)/'(&)+ / o{x)f{x)dx ( 23 ) 



1 f x 

P( x ) = g ~ i x J> a ( x ) = J o P^) dL 

Since a(x) < 1/8, it is easy to see that 

£ = ^r 4 + ^c(3)&- 3 + o(r 2 ), 

where ((x) is the Riemann zeta function. Therefore, b — > oo. As is known, ^(4) = 7r 4 /90. 

Similarly, in the right-hand side of pip , we can pass to integrals for sufficiently large 
values of I. 

In view of fTTTl) . the term 0(b 3 ) is greater than E\/N IniV and its calculation makes 
no sense. 

Similar estimates for N, as can be seen from Theorem 1 from [5J, also make possible 
the computation of the asymptotic (in b) terms 0(& -3 ) and 0(b~ 2 ) only. The leading 
(first) term of the asymptotics will be called the thermodynamical limit and the second 
one, the quasithermodynamical limit. 

Now let us consider a system of iV three-dimensional noninteracting oscillators of the 
same frequency ojq\ 

h 2 

- —A^ n (x) - Lu 2 \x\ 2 ^ n (x) = n* n (x), x E M 3 . (24) 

In order to obtain the leading term of the Stefan-Boltzmann law for this system of 
oscillators, we must set 

h = W 

here fi is the Planck constant, k is the Boltzmann constant, T the temperature, ujq the 
frequency, which equals ujq = cj \JV \ where c is the speed of light, and V is the volume. 
We assume that the oscillators are completely isotropic, so that the frequency co>o is the 



fi 



same in all directions. Oscillations of frequency u greater than uq do not exist at the given 
temperature, while the frequencies u < ujq give a considerably lesser number of variants 
and may be neglected. 

As a result, for the correction to the Stefan-Boltzmann law, we obtain a quasither- 
modynamical term of the form 

where F is the free energy, § = 7r 2 k 4 /(Q0h 3 c 2 ) is the Stefan-Boltzmann constant, V is the 
volume, c, the speed of light, h is the Planck constant, k is the Boltzmann constant, and 
T is the temperature. 

It has always been observed that the experimental values of the Stefan-Boltzmann 
constant are larger than their theoretical values. Now the reason for this discrepancy is 
clear: the correction term specified above explains this discrepancy. 

However, the main effect of quasithermodynamics occurs when there is no thermody- 
namocal phase transfer, but there is a quasithermodynamical one. It is precisely such an 
effect that was obtained by the author in the study of a Fermi-gas without the additional 
assumption on the existence of Cooper pairs. It was assumed that the number of particles 
N tends to oo, the particle interaction is pairwise, as in Helium 4, i.e, repulsive at short 
distances and attractive at large ones. At the same time, the potential V(ri — rj) is mainly 
a short distance one. The corresponding spectrum has the form 

E < = ^ + no - no), (26) 

where V(l) is the Fourrier transform of V(z). 

As N tends to infinity and the volume V tends to infinity in such a way that N/V — ► 
const, the potential tends to the 5-function and therefore tends to the spectrum of an 
ideal gas. Superfluidity arises only in the quasithermodynamical limit. In this limit, as 
a rule, the Landau curve arises. Thus, the phase transfer to the superfluid state occurs 
in the quasithermodynamics of this "model" . Note that this model does not contain any 
additional physical interactions: it is a "model without a model" , since the antisymmetric 
solution of the iV-particle Schrodinger equation with the ordinary pairwise interaction 
cannot be regarded as a "model" . This is an ordinary mathematical problem. 

Let us consider a Bose-gas. Bogolyubov proposed the following spectrum in this case: 



^=m/(|£ + V(Q) -nO 2 , (27) 

which also yields the Landau curve in quasithermodynamics. But unlike (1261) . in the 
thermodynamical limit it gives 



lun E, J[ — + V(0)) HO) 2 . 

jv/v-vconst, Af-^oo y y 2m J 

and superfluidity is preserved (the critical Landau speed in the thermodynamical limit 
is not zero, as it is in the case of the spectrum (126]) ). However, the photon part of the 
spectrum disappears in the thermodynamical limit and appears only in quasithermody- 
namics. 
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Remark. Bogolyubov's work may be rigorously founded under certain additional con- 
ditions, provided one considers the problem, as Bogolyubov did, on the three-dimensional 
torus. However, the passage to the limit from the torus to three-dimensional space is 
erroneous: an everywhere dense point spectrum appears in the limit. 

The author obtained a modification of formula (|27|) for the case in which the liquid 
flows through a capillary of radius r. The spectrum in that case has the form 



Ei 



l U^-kl)^ l - k j- V ^ 



v(i-k 2 )-v(k 2 y 



v(i + k 2 ) + v(k 2 )\ 2 fv(i-k 2 )+v(i+k 2 y 2 



1 1/2 



+ 2 ( a ( l i + {2 ~ 2k D + V(l - k 2 ) - V(k 2 )) 



x Ja 2 (l 2 - P f + 2(V(l + k 2 ) + v(k 2 )y 



(28) 



where 



h 2 

a=— , h = l + 2k 2 , 
2m 



while k 2 = 27m/r, where n = 1, 2, . . . . 

When r = oo, and hence k 2 = 0, we obtain the Bogolyubov formula (1271) . In the limit 
as r — > 0, k 2 — > oo, / = (Z , — k 2 ) (/ is the component along the direction of flow), we 
obtain 



UmE l0 = ^a^ + a\V(l )\l 2 . 



These results, just as the Bogolyubov formula, are valid in the thermodynamical limit 
with the quasithermodynamical correction. 
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